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Abstract: The dependence of the differential cross section da/dp± of inclusive heavy 
quark production in pp and pp collisions on the renormalization and factorization scales 
is investigated. The implications of our results for experiments at TEVATRON and LHC 
are discussed. In particular, it is shown that the NLO QCD predictions for tt production 
at the LHC based on the Principle of Minimal Sensitivity are by 30 — 50% higher than 
the standard ones. 
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1. Introduction 



Heavy quark production in hard collisions of hadrons, leptons and photons has been 
considered as a clean test of perturbative QCD. It has therefore come as a surprise that 
the first data on the bb production in pp collisions at the TEVATRON 7p collisions at 
HERA and 77 collisions at LEP2 have turned out to lie significantly and systematically 
above theoretical calculations. The disagreement between data and theory is particularly 
puzzling for the collisions of two quasireal photons at LEP2. 

The NLO QCD calculations of heavy quark production in these collisions depend 
on a number of inputs: a s , parton distribution functions (PDF) of colliding hadrons or 
(resolved) photons, fragmentation functions, masses of heavy quarks, and last but not 
least the choice of renormalization (RS) and factorization (FS) scales /i and M. 

The discrepancy between the DO [1] and CDF data [2] on bb production and theoretical 
calculations [3] has led to suggestions that it might represent the manifestation of effects 
of unintegrated gluon distribution function within the /^-factorization approach [4], or 
even a signal of supersymmetry [5]. In both cases the agreement between the TEVATRON 
data and the corresponding calculations [4, 5] is quite impressive. On the other hand, 
it has been argued in [6] that proper parameterization of the b-quark fragmentation 
functions also improves the agreement between CDF data and QCD calculations. There 
have furthermore been attempts to go beyond fixed order calculations by resumming the 
effects of large logs of the type ln(Sy4m 2 ) [7], where S denotes the center-of-mass energy 
squared, or ln(p^/m b ) [8], as well as the effects of soft gluon resummation important in 
the threshold region [9]. 

In [10] the sensitivity of NLO QCD calculations of the total cross section at t{PP —> bb) 
to the variation of the renormalization and factorization scales \x and M was investigated 
in some detail. It was shown that in order to arrive at locally stable results these two scales 
must be kept independent. Moreover, it turned out that in the TEVATRON energy range 
the position of the saddle point of the cross section af^°{pp — > bb; S, M,fi), considered 
as a function of \x and M, lies far away from the diagonal fx — M. Using the NLO 
prediction at the saddle point instead of the conventional choice /x = M = enhances 
the theoretical prediction in the TEVATRON energy range by a factor of about 2, which 
may help in explaining part of the excess of data over NLO QCD predictions. 

In this paper the analysis of [10] is extended to the differential cross section da/dp± 
of inclusive heavy quark production at the NLO. Bottom as well as top quark production 
will be considered in circumstances relevant for experiments at TEVATRON and LHC. 

2. Basic facts and formulae 

The basic quantity of perturbative QCD calculations, the renormalized color coupling 
ct s (fj), depends on the renormalization scale /i in a way governed by the equation 
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where for rif massless quarks (3q = 11 — 2n//3 and f3\ = 102 — 38n//3. The solutions 
of (2.1) depend beside fi also on the renormalization scheme (RS). At the NLO (i.e. 
taking into account first two terms in (2.1)) this RS can be specified, for instance, via the 
parameter A RS , corresponding to the renormalization scale for which a s diverges. The 
coupling «(/i) is then given as the solution of the equation 



& In ( -L- + C ln „ (2.2) 



where c = (47r/3o). At the NLO the coupling a s is thus a function of the ratio /x/Ars 
and the variation of the RS for fixed scale \x is therefore equivalent to the variation of [i for 
fixed RS. To vary both the renormalization scale and scheme is legitimate, but redundant. 
Let us emphasize that the choice of the RS is as important as the choice of renormalization 
scale. As setting the renormalization scale equal to some "natural physical scale" gives 
different results in different RS, the existence of such a "natural scale" is actually of no 
help in arriving at unique result. If not stated otherwise, we shall work in the conventional 
MS RS and vary the renormalization scale \i only. 

For hadrons the factorization scale dependence of PDF is determined by the system 
of evolution equations for quark singlet, nonsinglet and gluon distribution functions 



dS(M) 
dlnM 2 
dG(M) 
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where 



E(x,M) =^( ?i (x,M)+^(x,M)), (2.6) 
i=i 

g NS ,(x,M) = (q t (x,M) + q i (x,M)) - — E(x,M), Vi. (2.7) 

n f 

The splitting functions admit expansion in powers of a s (M) 

P tj(x , M) = ^Iplf ( X ) + (^l) 2 P$\ X ) + . . . , (2.8) 

where P^\x) are unique, whereas all higher order splitting functions Pjl{\j > 1 depend 
on the choice of the factorization scheme (FS). The equations (2.3-2.5) can be recast 
into evolution equations for qi(x, M), q^x, M) and G(x, M). In all calculations presented 
below we took rif = 4 in the case of bottom quark and rif = 5 for top quark. The 
RG equation (2.2) was solved numerically in order to guarantee correct behaviour of its 
solution for small renormalization scales. 
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3. General form of the cross section 



According to the factorization theorem, the differential cross section for inclusive QQ 
production in hadronic collisions has the form (suppressing the specification of the de- 
pendence on rriQ and y/S) 

da(p±,y) \- [f A , . , da^x^M.p^y) 
dydp± = J J teitetfMxu M)f j/B (x 2 , M) dydp± ' (3 ' 1} 

where y and pj_ denote the rapidity and transverse momentum of a heavy quark Q, A 
and B stand for the beam particles and the partonic hard scattering cross sections admit 
expansions in powers of a s (fj) 



d0-ij(xiX2,M,pj.,y) _ 2 (2 )( ^ j 3 ^ ( 3 ) 



= al(fi)c^\xix 2 ,p±,y) + a 3 s (fx)c\f (x!X 2 , M, fi,p ± ,y) + ■ ■ ■ . (3.2) 



Note that whereas the factorization scale scale M appears as argument of both the PDF 
fi/A,fj/B and parton level hard scattering cross sections a^, the renormalization scale \x 
enters only if we calculate the latter as perturbative expansion in powers of a s . Both 
the renormalization and factorization scales are auxiliary, unphysical quantities, on which 
observables like the cross section (3.1) cannot depend, provided both the expansions (3.2) 
and (2.8) are taken to all orders. In such a case any choice of /x, M gives the same result for 
(3.1). If, however, these expansions are considered to a finite order, as is in practice always 
the case, the expression on the right hand side of (3.1) develops nontrivial dependence 
on these two scales and their choice therefore matters. Moreover, as the renormalization 
scale reflects ambiguity in the treatment of short distances, whereas the factorization scale 
comes from similar ambiguity concerning large distances, these two scales should be kept 
as independent free parameters of any finite order perturbative approximations [11]. 



4. Choice of scales 

For a given p±, the result of the evaluation of (3.1) to the NLO can be represented graph- 
ically in the form of two-dimensional surface over the plane /i, M. In this representation 
any of the points on such a surface corresponds to one, in principle equally legitimate, 
choice of these two scales. By constructing this surface we get quantitative idea of the 
stability of our calculations. Specifically we can look for the presence of saddle points, 
where the NLO results exhibit locally the independence of /i and M, which is a necessary 
feature of all-order calculations globally. The choice of scales based on this criterion, 
called Principle of Minimal Sensitivity (PMS) [12], is one of the two general approaches 
to scale setting, the other being the Effective Charges method [13]. This latter approach 
selects those points M E for which the LO and NLO approximations of (3.1) coincide 
and the perturbative expansion thus has the "fastest apparent convergence" . Graphically 
this choice of scales is given by the intersection of LO and NLO surfaces representing the 
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cross section da/dp±(fi, M). As shown in [10] for the total cross section the resulting 
points form a curve which passes very close to the saddle point, if such a point does exist. 
As will be shown below the same holds for differential cross section da/dp^. 

In this connection we want to recall the simple but often overlooked fact that the 
conventional procedure of setting \i = M and identifying it with some "natural physical 



scale" (in our case usually mj_ = yrrtQ + p 2 ± ) does not actually resolve the underlying 

ambiguity, because it implicitly assumes that we are working in the MS RS. There is, 
however, nothing "natural" on this (or any other) RS, as the MS is just one of the infinite 
number of equally legitimate choices of the RS. As emphasized in Section 2, working in 
a different RS of the couplant a s (fi), the same choice of "natural" renormalization scale 
\i leads to a different result for the cross section (3.1). This fact also implies that it is 
meaningless to estimate the "theoretical uncertainty" related to freedom in the choice 
of the renormalization scale /z by varying it in some interval around the "natural" scale 
because the same procedure leads in different RS to different bands of theoretical predic- 
tions. The standard choice of MS RS is often justified by noting that it leads to small 
coefficients in perturbative expansions of physical cross sections. But if the apparent rate 
of convergence is the basic criterion for the "proper" choice of the scales than the points 
preferred by the EC approach, rather than those defining the MS RS should be selected. 

On the other hand, although the positions of the saddle points do also depend on the 
chosen RS, the values of the cross section (3.1) at these saddle points do not! Contrary 
to the scale setting methods relying on the existence of "natural physical scale" the PMS 
approach therefore results in a unique prediction. For predictions based on the method 
of effective charges [13] the same uniqueness holds for any fixed factorization scale M 
as well. In the case both /i and M are treated as independent free parameter the EC 
approach does not yield a unique prediction but the whole continuum of values reflecting 
the fact that the intersection of LO and NLO surfaces defines a curve 1 . 

However, independently of the strategy one prefers for choosing the renormalization 
and factorization scales, it is undoubtedly useful to have quantitative idea of the sensitivity 
of our calculations to the variation of these unphysical parameters. 

5. Results 

The cross section (3.1) was computed to the NLO using the code of the ref. [3] together 
with CTEQ6.1 set of PDF (with A4 = 326 MeV). As no restrictions on the rapidity y 
of the produced heavy quark was imposed in this study, the integration over the whole 
kinematically allowed region of y was performed. To see how much our conclusions depend 
on the knowledge of PDF, some of the calculations were done also for the GRV94 set. The 
difference between the CTEQ6.1 and GRV94 results, should, however, not be taken as 

1 This statement should be taken with some reservation as there is in fact no guarantee that the saddle 
points and/or the mentioned intersection do, in a given kinematic region, exist. 
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an estimate of their dependence of the chosen PDF. We have deliberately taken the now 
obsolete GRV94 set in order to illustrate the independence of the basic features of the 
renormalization and factorization scale dependence of do" NLO /dpj_ on the selected PDF. 

The numerical integration used in [3] employs standard Vegas method. For the eval- 
uation of the LO contribution 4 iterations each with 4000 points were used, whereas for 
the NLO term these numbers were increased to 6 iterations each with 8000 points. 

5.1 bb production 

For inclusive bb production the differential cross section da/dp± was computed for pp 
collisions at = 64, 630 and 1800 GeV and for pp collisions at v^S* = 14 TeV. For 
pp collisions the first energy does not correspond to any real experiment, but we have 
included it to make contact with the earlier studies of this kind [14]. In all calculations 
we set nib = 4.95 GeV. 

We start with pp collisions at = 64 GeV. In Fig. 1 the quantity (3.1), integrated 
over the rapidity y, is displayed as a function of \i and M for several values of p±. On each 




Figure 1: (Left) Positions of saddle points, b-quark cross section and ratio (5.1) for CTEQ6.1 
(green) and GRV94 (black) PDF at y/S = 64 GeV. Red line shows conventional choice fj, = M = 
\Jp\ + mf- (Right) Surfaces plots for some four typical values of p±, calculated with CTEQ6.1 
PDF. Yellow circles mark saddle points and yellow triangles conventional choice. 
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surface the point corresponding to the standard choice [i — M — \Jp\ + rn% is marked 



by a yellow triangle. Each of the surfaces in Fig. 1 exhibits a saddle point, marked by 
a yellow circle. The p± dependence of the coordinates (/* s (p_l), M s (p±)) of these saddle 
points, displayed in the upper left part of Fig. 1, shows the first indication of the disparity 
between the renormalization and factorization scales: saddle points occur at locations 
clearly away from the diagonal /j, — M and always for ji s < M s . Moreover, for all values 
of p± these saddle points lie almost exactly on the intersection (not shown) between 
the LO and NLO surfaces, which defines the points preferred by the Effective Charge 
Approach. We shall plot the corresponding contour plot, illustrating this claim, for the 
phenomenologically most interesting case of tt production at LHC. In the middle left part 
of Fig. 1 the distribution da/dp± computed for the standard choice [i — M— \Jp\ + 
is compared to that corresponding to the saddle points. To quantify their difference, we 
plot in the lower left part of Fig. 1 the ratio (setting mq = m b ) 



The difference between the results obtained with CTEQ6.1 and GRV94 PDFs is negligible 
as far as the position of saddle points is concerned, but CTEQ6.1 gives, for both the PMS 
and the conventional choice of scales, slightly higher values of the cross section. In the 
region p± < 10 GeV saddle-based NLO results are by a factor of about 1.5 higher than 
the conventional ones. 

As we go to higher energies the relation between the conventional and saddle-based 
results changes. At V~S = 630 GeV (see Fig. 2) we find that 

• In the interval p± < 100 GeV the saddle point lies far away from the diagonal 
li = M. 

• For p± < 30 GeV the renormalization scale at the saddle point /i s ~ 1 GeV, dipping 
to very small value around 0.6 GeV for p±_ = 16 GeV, where the ratio (5.1) peaks 
at about 10. 

• Stability of the NLO predictions improves for p ± > 100 GeV, where the saddle point 
(/i s , M s ) moves closer to the conventional choice /i = M = >Jml +p 2 ± and the ratio 
(5.1) approaches unity. 

• The shapes of the conventional and saddle-based p±_ distributions differ substantially, 
particularly in the region 10 < p± < 40 GeV. 

At this energy and in the experimentally accessible region p± < 30 GeV the typical parton 
fractions come from the region x > 0.03 and thus neither the effects of low x and/or large 
p_i_/m b are likely to change the fixed order predictions significantly. 

For pp collisions at = 1800 GeV (see Fig. 3) the saddle point disappears for 
P± ^ 10 GeV and reappears only for p ± > 80 GeV. The value of the renormalization 



r(p±) 



da(/i s , M s ,p ± )/dp ± 



(5.1) 
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Figure 2: The same as in Fig. 1 but for = 630 GeV. 



scale at the saddle point, fj, s , remains small, dropping below 1 GeV as p± approaches 10 
GeV. The saddle point remains close to the intersection of the LO and NLO surfaces even 
as the saddle point approaches the unphysical region, though in this case the distance 
between the saddle point and the mentioned intersection increases. For p± < 5 GeV 
the ratio (5.1) stays constant at about 2, but then rapidly rises, diverging at p± = 10 
GeV. Only for very large and experimentally inaccessible region p± > 160 GeV do the 
conventional results approach the saddle-based ones. At these large transverse momenta 
the resummation of the large logs of the type \n(p±/rrib) must certainly be performed 
to improve the NLO calculations. However, in most of the p± range accessible at the 
TEVATRON, i.e. p± < 50, this effect is likely to be small. For moderate p± the NLO 
QCD predictions of inclusive bottom production are thus highly unstable. For very small 
roughly p± < 10 GeV, one may argue that because the relevant parton fractions are 
typically 0.01, the approach [4] using the unintegrated PDF within the /c^-factorization 
may be superior to the present fixed NLO results. It should, however, be kept in mind that 
the agreement with TEVATRON data relies on the use of a particular set of unintegrated 
PDF (which are still not well-constrained from data) and the appropriate choice of the 
factorization scale. In the case of /c^-factorization there is, however, no mechanism for the 
cancelation of the factorization scale dependence between the unintegrated PDF (which 
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Figure 3: The same as in Fig. 1 but for VS = 1800 GeV. 



depend on it) and the off-mass-shell hard scattering cross section as the latter is taken to 
the LO only and is thus independent of M. 

For the proton-proton collisions at the LHC the interval where the NLO QCD pre- 
dictions exhibit no stability region and behave as a monotonously decreasing functions of 
the renormalization scale \x extends between 10 < pj_ < 900 GeV. In the lower part of 
this range the partons (mostly gluons) involved in hard collision have momentum frac- 
tions typically of the order of 10~ 3 , for which the calculations [4] based on unintegrated 
PDF is likely to be more appropriate. On the other hand, for large transverse momenta 
the resummation of logs Inp^/m^ is definitely necessary to arrive at a reasonably stable 
result. Compared to these effects the threshold ones are likely to be negligible. 

5.2 tt production 

We have computed da/dp± for pp collisions at y/S = 1800 GeV and for pp collisions at 
\/S = 14 TeV, setting in both cases m t = 175 GeV. 

For the former case the results, presented in Fig. 4, exhibit saddle points for p± < 
200 GeV. As these saddles are broad, the points corresponding to the standard choice 
H = M = m t are not far away and the saddle-based and standard results are close, typ- 
ically to within less than 10 %. Moreover, contrary to the case of bb production, the 
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saddle points are located at large renormalization scale fi s > 80 GeV. At the TEVA- 




Figure 4: The same as in Fig. 1 but for tt production at \fS = 1.8 TeV. 



TRON about ~ 85% of the cross section comes from light quark qq annihilation with the 
remaining roughly 15% from gluon-gluon fusion. The typical parton fractions x q , Xq ^ 0.2 
come from the region where quark distribution functions are very well (with the precision 
better than about 3%) constrained by existing data, but the gluon distribution function 
suffers from much larger uncertainties. Recent error analysis by both CTEQ and MRST 
groups 2 show that just above xq — 0.2 this uncertainty starts to increase quite rapidly 
reaching 30% above xq = 0.3 for CTEQ or xg = 0.4 for MRST, respectively. Despite 
the fact that gluon-gluon channel contributes just a small part of the cross section, the 
overall uncertainty may thus be comparable to or even larger than that resulting from 
the scale ambiguity. This fact was noted in [15] 3 which contains detailed investigation 
of the dependence of NLO QCD predictions, including the soft gluon resummation, on 
the variation of the scales as well as on choice of parton distribution functions. The soft 
gluon resummation reflects the fact that at TEVATRON the tt production threshold is 

2 The numbers quoted below concerning the accuracy of PDF were obtained with the help of Durham 
HEPDATA PDF interactive facility [17] using CTEQ6 and MRST E-sets. 

3 Where the theoretical uncertainty was defined in a standard manner by changing the common scale 
/i = M by a factor of 2 around m±. 
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not far below y/S. Unfortunately, as the calculations of [9], used in [15], are available for 
the total it cross section only, we cannot investigate their implications for the transverse 
momentum distributions considered in the present paper. Taking as an estimate of their 
importance the number given in [9] for the total cross section cr tot (pp — > tt), i.e. the 
enhancement of the standard NLO results with [i = M = m t by about 4%, we expect 
them to be comparable to the difference between the standard and saddle-based NLO 
results. However, as the scale dependence (assuming \i = M) of the resummed results 
for a(pp — > tt) is reduced substantially with respect to standard NLO calculations [9], it 
would be quite interesting to see the effects of the soft gluon resummation on the surfaces 
in Fig. 4. 




Figure 5: The same as in Fig. 1 but for tt production in pp collisions at \/~S = 14 TeV. 



Though non-negligible, all the uncertainties of higher order QCD corrections discussed 
above are nevertheless substantially smaller than the current experimental errors even for 
the total cross section. It is difficult to quote just one number as the results from various 
channels and from DO and CDF experiments differ sizably (see Figs. 5 and 7 of a recent 
review [16]), but a conservative estimate of the experimental uncertainty for the total cross 
section is 30 — 50%, depending on the way one adds statistical and systematical errors 
(the number in [16] reads a to t(pp tt) = 8.2t2^(stat)~^{'^(syst)). The only conclusion 
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Figure 6: Left: the standard way of estimating the uncertainty of NLO calculation (lower solid 
black curve) by evaluating the cross section in the band between \i = M = mr/2 (upper dotted 
curve) and \i = M = 2mr (lower dotted curve), where mj- = y^mf The saddle-based 

results are represented by the upper solid green curve. Right: contour plot corresponding to 
the surface plot in Fig. 5 where the dotted curves show the choice of scales preferred by the 
approach of [13]. 



one can thus draw from the analysis of the TEVATRON data on a tot (pp — > tt) is that 
they are consistent with QCD calculations, with or without the soft gluon resummation. 

At the LHC and even in the initial low luminosity period, tt pairs will be produced 
abundantly, with millions of such pairs produced per year by each of the two main exper- 
iments. The gluon-gluon channel contributes most of the cross section, except for very 
large p± > 700 GeV, there the light quark-gluon channel becomes of comparable impor- 
tance. The qq channel contributes less than 10% throughout the experimentally accessible 
regions. The precision of the measurements of tt production will thus be limited primarily 
by the uncertainty of the luminosity, which is expected to be in the range 5 — 10% [18]. 

As shown in Fig. 5, for tt production at the LHC the saddle points exist for all 
accessible values of p±, but they are located further away from the points corresponding 
to the standard choice, primarily due to the fact that the renormalization scale at the 
saddle points, /x s , remains below 100 GeV. This by itself is not so important, what matters 
is the fact that the difference between the saddle-based results and the conventional is 
quite sizable. For p± < 100 GeV and the CTEQ6.1 set of PDF the former are higher by 
a factor of 1.3 — 1.5. This is by about a factor of two more than the upper limit of the 
theoretical band based on the standard way of estimating the "theoretical uncertainty" of 
NLO calculations, shown by the lower (black) solid curve in the left part of Fig. 6. This 
enhancement is phenomenologically relevant since at LHC the typical parton fractions 
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will be by a factor of seven smaller than at the TEVATRON and thus in the region where 
both light quark and gluon distributions functions are known with accuracy better than 
about 3% [17]. The sensitivity of NLO QCD predictions to the choice of PDF is thus 
expected to be substantially smaller than the difference between the standard and saddle- 
based NLO calculations, with or without the soft gluon resummation. The effects of the 
resummation of large logs coming from the region of small x and/or large p±/m t are also 
unlikely to be important in the experimentally accessible region. 

As demonstrated by the four contour plots in Fig. 6 the saddle points are located 
close to the curves representing the intersections of the LO and NLO predictions, the 
feature mentioned already in the case of bb production at v^S* = 64 GeV. Note moreover, 
that for all p± these intersections follow closely the ridge on both sides of the saddles 
implying that out of all the points preferred by the Effective Charges Approach [13] the 
one closest to the saddle gives the lowest value of the cross section. 

6. Summary and conclusions 

We have investigated the renormalization and factorization scale dependence of inclusive 
bb and tt production in pp and pp collisions as a function of p±. For bottom quark at 
64 GeV and top quark at TEVATRON and LHC energies, the cross sections were found to 
exhibit regions of local stability around the saddle points for all experimentally relevant 
values of p±. In all these cases the standard choice of scales /j, — M — \Jp\ + m 2 yields 
results which are lower than those at the saddle points. The difference is smallest for tt 
production at TEVATRON energie, where it amounts to about 10%. For pp collisions at 
the LHC the enhancement of the saddle-based predictions over the standard one increases 
to to 1.3 — 1.5, depending on p±. As at LHC top quark production cross sections will 
be measured with a precision of about 5 — 10%, and hopefully even better, theoretical 
predictions should be of comparable accuracy. In view of this fact the latter enhancement 
is phenomenologically relevant. This is our most important conclusion. 

For smaller ratio m/V~S, characteristic for the bottom production at TEVATRON 
and LHC, the saddle points are located far away from the diagonal [i = M and the corre- 
sponding saddle-based results are substantially higher that the standard ones. For some 
regions of p± there is no region of local stability at all. For small values of p± the process 
is dominated by very small partonic fractions and thus the use of unintegrated PDF may 
be more appropriate. For large p± > 50 GeV, on the other hand, the resummation of the 
large logs of the type \npj_/rrib should certainly precede any scale analysis. 
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